Theory of neutral nitrogen- vacancy center in diamond and its qubit application 
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The negatively charged nitrogen- vacancy defect (NV)^ in diamond has attracted much attention in recent 
years in qubit and biological applications. The negative charge is donated from nearby nitrogen donors that 
could limit or stem the successful application of (NV)~. In this Letter, we unambiguously identify the neutral 
nitrogen- vacancy defect (NV") by ab initio supercell calculations. Our analysis shows that i) the spin state can 
be selectively occupied optically, ii) the electron spin state can be manipulated by time-varying magnetic field, 
and iii) the spin state may be read out optically. Based on this NV" is a new hope for realizing qubit in diamond 
without the need of nitrogen donors. 



Realization of qubits is of extremely high importance be- 
cause that can be used in quantum cryptography, quantum 
optics and quantum computing. One of the most promis- 
ing candidate is the negatively charged nitrogen-vacancy de- 
fect (NV~) ill \M that can operate at room temperature 
iS B II S Sli HE El- The negative charge is donated 



from nitrogen substitutional (Ns) 11121 1131] . It has been very 



recently shown that the major source of the decoherence of 
the electron spin of the negatively charged NV center is the 
electron spin bath of nitrogen substitutionals in type-/6 dia- 
mond, that can be eliminated only at low temperature using 
a giant magnetic field |14]. This is a serious limiting factor 
towards the development for practical applications. In addi- 
tion, it has been recently demonstrated that nanometer-sized 
diamond particles containing NV^ are useful as fluorescent 
biomarkers for in vitro imaging applications 115, T^. How- 
ever, the requirement of the extra charge on NV defect from 
Ns could be critical in biomarker applications, where the size 
of the nanodiamond is reduced to some nanometers in diam- 
eter 11711 . One solution of this problem would be to apply 
the paramagnetic neutral NV defect (NV°) which produces 
strong photoluminescence at 2.156 eV (similar to NV^) and 
it does not require an extra electron. The NV" center is well 
identified Ill2l1 while its strong connection to its neutral coun- 
terpart is well established experimentally |18, 23] ■ How- 



ever, the overall knowledge about NV° is scarce. Very re- 
cently, an = I electron paramagnetic resonance (EPR) cen- 
ter has been found in photo-excited diamond doped by ^^N 
1 2111 . A sizeable ^^N hyperfine constants have been detected 
in the EPR measurements and it was proposed that the signal 
was originated from one of the excited states of NV" Kill . We 
emphasize that the identification of the EPR signal of NV" is 
a key step in order to apply it for qubit, and to trace it magnet- 
ically in biomarker applications. 

In this Letter we i) unambiguously identify the EPR signal 
of NV" ii) provide a detailed spin density distribution around 
the defect; the results indicate that the defect states are well- 
localized and the electron spin can be decoupled from the 
spin bath of ^^C nuclei iii) analyze the electronic structure 
by group theory explaining the photo-excitation of NV°; this 
shows that the A/s = ±1/2 sublevels are selectively occupied 
in the photo-ionization process and appropriate microwave 
magnetic field can be applied to choose also selectively ei- 



ther A/s = +3/2 or A/s = —3/2 which finally can be used 
as qubit iv) we propose that the A/s = +3/2 or A/s = —3/2 
states may be read out optically during the emission process 
from the excited state to the ground state. This gives a hope 
for realization of qubit by NV'^ defect. 

We employed density functional theory calculations with 
local density approximation (LDA) using a large, 512-atom 
simple cubic diamond supercell. We used two different codes: 
the geometry of the defect was optimized with the VASP 
code 1 22] while the hyperfine tensors of NV° were calculated 
by the CPPAW code {2% Both VASP and CPPAW codes 
apply the all-electron PAW method and plane wave basis set. 
We used a cut-off of 30 Ry and F-point for k-point sampling. 
Other details and references can be found in our previous pub- 
lication 1(24}] where we could successfully describe the nega- 
tively charged NV center and can reproduce weU the experi- 
mental hyperfine data. 

The NV defect has €3^ symmetry, if no reconstruction oc- 
curs. We use defect-molecule and group theory analysis below 
as for NV" (c.f. [24] and references therein). Four dangling 
bonds point to the vacant site resulting in ai (l)^ai (2)^e^ con- 
figuration for the neutral defect, where ai(2)^ and the dou- 
ble degenerate appear in the fundamental band gap. In 
this configuration, the system is basically Jahn-Teller unstable 
because the degenerate state is only partially filled by elec- 
trons. The static Jahn-Teller effect would result in Cih sym- 
metry of the defect. PL studies indicate that the system has 
Cat, symmetry {^E state) with possibly exhibiting dynamic 
Jahn-Teller effect that inhibits the EPR detection of the "^E 
ground state i25ll . The Cih symmetry configuration and the 
"^E ground state of €3^ symmetry can be described by single 
Slater-determinants, so we could address this issue directly 
at LDA level. The lowest excitation can be obtained by pro- 
moting one electron from ai (2) level to e level resulting in 
ai(l)^ai(2)^e-^ configuration. The possible excited states of 
ai(l)^ai(2)^e^ configuration have "^Ai, and "^E multi- 
plets. The A/s = | state of multiplet can be described 
by a single Slater-determinant by simply aligning all the elec- 
trons spin-up on the ai(2) and e levels, thus the spin density 
can be determined by usual LDA calculation. We note that no 
Jahn-Teller effect occurs for * state, so we considered only 
Cat, symmetry for this state in the calculations. We focused 
our research on the ground state of Ciu and Cat, symmetries. 
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and the excited state that are relevant in the recent EPR 
study ||21|1 . We note here that the nitrogen dangling bond is 
hybridized into Oi (2) but not into the e defect states, there- 
fore, negligible spin density is expected for the '^E ground 
state but considerable spin polarization may be expected for 
the excited state. 

First, we investigated the ground state of NV*^. Since Jahn- 
Teller effect can occur, we distorted the symmetry to Ci^, and 
allowed the atoms to relax to find the energy minimum within 
LDA. We found that the system conserves its Cui symmetry 
and it does not form the C^y configuration. The single oc- 
cupied e-level in the gap is split by 0.3 eV resulting in an 
occupied a' state and an unoccupied a" state of the spin-up 
electrons. The neighbor N and C atoms relaxed outward from 
the vacant site. The nitrogen atom remained very close to the 
(111) axis even in Ci/i symmetry, while one of the carbon 
atoms is considerably closer to the vacant site than the other 
two. Next, we constrained the system to preserve the C^y 
symmetry during the geometry optimization in order to cal- 
culate the total energy of state. We found that it is about 
0.09 eV higher in energy than the Jahn-Teller distorted one. 
Then, we allowed the system to relax without symmetry con- 
straints starting from the optimum C^y geometry. The atoms 
automatically relaxed to the Ci^ symmetry. This result indi- 
cates that symmetry is a local maximum. This finding 
seems to contradict the PL spectrum which shows ^E ground 
state within C^y symmetry. However, we carried out quasi- 
static calculations at K, and we did not take the vibronic 
states and temperature effects into account. This result shows 
clearly a multi-valley potential surface for the ground state of 
this defect as was already hinted by Davies |25(1 . Our calcu- 
lated Jahn-Teller energy (0.09 eV) is very close to the esti- 
mation of Davies (0.14±0.07 eV) based on the luminescence 
measurements |25]. There are three equivalent Cih configura- 
tions around the vacant site rotated by 120 degrees about the 
Cat, -axis with each possessing the global minium. At LDA- 
level, we can estimate the upper limit of the energy barrier 
between the three global minima, that is «0.09 eV. There may 
exist more favorable path between these global minima than 
through the Cs^ configuration, so the actual barrier energy 
may be even lower than that. Due to the strong C-C bonds 
in diamond <0.09 eV energy difference can be gained even at 
very low temperature because this is about the zero point en- 
ergy of the phonons with the highest energy (<0.083 eV). This 
results in a motional average of the single electron between the 
three equivalent Ci^ configurations, showing an effective C^y 
symmetry. It was speculated 12 ll 12511 that the dynamic Jahn- 
Teller effect is responsible for the missing EPR signal of the 
'^E ground state. Our calculations support this assumption. 

We also investigated the '*A2 excited state. This state has in- 
deed much higher energy by about 0.86 eV compared to that 
of the low-symmetry ground state. We note that the LDA to- 
tal energy differences should not be directly compared to the 
measured transition energies as previously discussed in Ref. 
[24. In this state the ai(2) defect level is polarized, thus the 
nitrogen atom is polarized in contrast to ^E state. This re- 



TABLE I: The calculated principal values of the hyperfine tensor 
(columns 2 to 4) compared to the known experimental data (columns 
5 to 7) in MHz. The experimental data on ^""N is taken from Ref.[2ll 
The ^"^C hyperfine signal could not be resolved in that experiment. 
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A22 
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.MP 


^33 


^=*C(3x) 


-23.4 
60.6 


-23.4 
61.0 


-39.0 
126.2 


-23.8(3) 


-23.8(3) 


-35.7(3) 



suits in considerable rearrangement of the atoms around the 
vacancy, namely, the nitrogen moves closer to the vacant site. 
We calculated the hyperfine tensors of atoms for the optimum 
geometry and compared to a recently found EPR center as ex- 
plained in the introduction 112 ill (see Table IB . The agreement 
between the calculated and measured hyperfine signal is ex- 
cellent. In addition, the calculated binding energy of the NV° 
complex (w3.0 eV) shows a high thermal stability, and its (- 
/O) occupation level is at about 2.0 eV above valence band 
edge, in line with the experimental observations [18]. Thus, 
we identify the EPR signal of NV°. We provide the hyperfine 
data of the three nearest carbon atom near the vacant site that 
has the largest hyperfine interaction. These carbon atoms may 
be measured by future EPR experiments when the signal to 
noise ratio can be reduced there. 

After identification of NV*^ defect we discuss its possible 
role in spin physics. NV^ was successfully used to realize 
qubits igi [Tol]. Second order correlation and EPR measure- 
ments were employed to detect individual NV^. The result- 
ing spin-echo signals show a rapidly oscillating function en- 
veloped by a more slowly oscillating function [9]. The authors 
of Ref. 9 proposed a theory to explain this signal, and they 
concluded that the fast modulation frequency is due to the ef- 
fective magnetization density of the electron spin felt by the 
^^C nucleus, which is the same as the hyperfine interaction. If 



too many C nuclei are involved in the process that will lead 
to fast decoherence of the spin-echo signal |9], therefore, the 
knowledge of the spin density distribution is crucial. The situ- 
ation is complex for the ground state as was discussed above. 
At the measurement temperature an effective C^v symmetry 
may be detected as the motional average of the three equiva- 
lent Cih configurations. It is difficult to handle this situation 
by quasi-static simulation. We simulate this average simply 
by taking the optimum ^E ground state within C31, symmetry. 
We found that the overall picture is very similar to what was 
found for NV^ [27]. In the ground state the spin density is 
mainly localized on the three nearest neighbor carbon atom of 
the vacant site, and the spin density decays fast as a function 
of the distance from the vacant site. From this point of view, 
a potential qubit of NV" would be weakly coupled to its en- 
vironment and, therefore, it might present the right coherence 
properties like NV~. 

However, the localized spin density is not the only require- 
ment to produce qubits. One needs to generate a superposi- 
tion state and read-out the qubit states. In the case of NV^ the 
A/s = of the triplet ground state can be optically pumped 
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which has much smaller hyperfine interaction with the prox- 
imal ^"^C nucleus than the Ms = 1 state [9]. This effect was 
responsible for the collapse and revival of the electron spin 
coherent state. The Ms states could be read-out also optically 
after the measurements by using the fact that they have differ- 
ent fluorescence rates |9]. 

In NV*' defect the optical pump can change the S = ^ state 
to 5* = I excited state, and by switching off the light exci- 
tation this can be transformed back to the S = \ state. The 
nature of the spin-flip process will be discussed shortly. The 
group theory analysis tells us that no spin-orbit coupling arises 
for the '^A2 excited state. However, the spin-spin interaction 
is active, which can be given by the following effective spin- 
Hamiltonian for this particular system: Hss — D'{S^ — 5/4), 
where D' is the zero-field constant and Sz has the Ms eigen- 
value. So, the 4x degenerate state will split to two dou- 
ble degenerate states due to spin-spin interaction resulting in 
the lower lying Ms = 1/2;— 1/2 and the upper lying Ms = 
3/2; —3/2 levels (see Fig. These levels are separated by 
2 X D' = D. D K, 1685 MHz was experimentally measured 
by EPR iiit]. Acording to Davies PL analysis Q ~ 2.2 eV 
excitation occurs between ground state to the "^Ai excited 
state. We propose that the Ms = 1/2; —1/2 electron states of 
the '^Ai excited state can relax to the Ms ~ 1/2; —1/2 sub- 
levels of ^A2 state with a finite probability instead of relaxing 
back to the ground state. While the spin-orbit interaction is not 
active for '^Ai state itself the axial spin-orbit interaction cou- 
ples the ^Ai states selectively with the Ms — 1/2; —1/2 
states. That is the source of the spin-flip process which may 
be further mediated by phonons in order to satisfy the energy 
conservation of this transition. Indeed, the treshold excitation 
energy of the EPR signal, 2.2(1) eV jUJ, is larger than the 
ZPL energy of 2. 156 eV, and it has the highest intensity using 
excitation energy of about 2.5 eV. The probability of the relax- 
ation process between the original and ^Ai (S = ^) states 
should be higher than for the original (S = |) state. One 



may estimate from the known data of the NV^ center |126 
that the duration of the direct transition process is about 10 ns, 
while for the spin flip process it could be about 30 ns. There- 
fore, the duration of the optical pumping can be relatively long 
in order to arrive at Ms = 1/2; —1/2 sublevels of '^^2 state 
from the "^E ground state. The typical EPR condition for ''A2 
state is shown in Fig.[T| The small external constant magnetic 
field (B) will spHt the Ms = 1/2; -1/2 sublevels lowering 
the energy of Ms = — 1/2 state. Finally, we predict that dur- 
ing the optical pumping the Ms = 1/2; —1/2 sublevels of 
'^A2 state will be selectively populated with somewhat higher 
probability for Ms = —1/2 because it is lower in energy. 
Then, the varying microwave magnetic field induce the tran- 
sitions AIs ~ ±1/2 ^ ±3/2 in the EPR measurements. This 
scenario can explain all the photo-EPR findings in Ref. |2 lJ. 

The question arises what kind of entity can be used as a 
qubit from NV'^ defect. The first choice it to use simply the 
5-state. By optical pumping the S-state of NV° defect can be 
transformed from 1/2 to 3/2, and by switching off the li ght i t 
can be transformed back from 3/2 to 1/2 with t ^ 1 s 112 III . 



However, it is not probable that coherent state can be achieved 
for the ground state as we showed that it exhibits a dynamic 
Jahn-Teller effect causing a rapidly varying effective magnetic 
field for the ^^C nuclei around the NV° defect. Furthermore, 
the time-averaged effective hyperfine interaction of both states 
are very similar. This is also disadvantageous. 

Another possibility is to use the Ms sublevels of the '^A2 
state as qubit. As explained above by optical pumping one can 
select the Ms = 1/2; -1/2 sublevels of "^As state of the NV" 
defect with almost equal probabihty. That is not applicable for 
qubits. However, one can selectively set either Ms = +3/2 
(with energy of hvi) or Ms = —3/2 {hi'2) states by applying 
a TT pulse H to induce the EPR transition (AAfs — ±1) (see 
Fig-[D- The Ms ~ ±3/2 states may be used as qubit. It is 
clear in the above mentioned scenario that i^i can only be asso- 
ciated with Ms = +3/2 while 1^2 only with Ms = -3/2. We 
would like to emphasize that these metastable Ms — ±3/2 
states are extremely long living (>1 /is) because i) the '^A2 
state is the only 5 = 3/2 state, so there is no way for ra- 
diative recombination for this state ii) the Ms = ±3/2 states 
are not coupled to the ^Ai excited state at all that hinders 
the recombination of these states to the "^E ground state via 
"^Ai excited state. This is unique compared to the NV^ center 
1.26.1 and suggests a longer lifetime of these metastable states 
than for the measured lifetime of the singlet metastable state 
in NV- center (300 ns) fl^. We note that the Ms = ±3/2 
states are very weakly coupled to the ^E ground state by non- 
axial spin-orbit interaction. However, the non-axial spin-orbit 
interaction is very small (as assumed for NV~ center |26l]) 
and it is too far in energy in order to mediate this process 
by phonons. So, the probability of this decay is very small 
ensuring the very long lifetime of these states. The coher- 
ent state between these Ms = ±3/2 states and the proximate 
^■^C nuclei can certainly be maintained similar to NV~ 
The readout process could be very simple: by applying a tt 
pulse again the Ms — ±3/2 states scatter to ''A2 Ms — ±1/2 
states that through the spin-orbit coupling can go back to the 
"^Ai excited state, finally by radiative recombination to the ^E 
ground state. Therefore, the spin qubit state can be readout op- 
tically. This process can be tuned by applying the appropriate 
constant magnetic field to split the levels and the microwave 
magnetic field (tt pulse) to induce transitions between the lev- 
els. 

In this work we investigated the neutral nitrogen-vacancy 
defect in diamond in detail by ab initio LDA supercell calcu- 
lations. We showed that the defect indeed shows the dynamic 
Jahn-Teller effect for the ground state. We identified the re- 
cently found EPR center |21], as the '^A2 excited state of the 
neutral nitrogen-vacancy defect. That EPR center can be used 
to trace the NV"-contained nanodiamonds magnetically. We 
found that NV" is a promising candidate for realizing qubits 
in diamond without the need of nitrogen donors. 

AG acknowledges support from Hungarian OTKA No. K- 
67886. The fruitful discussion with Jeronimo Maze is appre- 
ciated. 
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FIG. 1: (Color online) The process of the manipulation of a 
qubit in NV" defect. Straight red arrow: radiative recombina- 
tion. Curved green arrow: spin-orbit coupling (Hso) possibly me- 
diated by phonons. Grey dotted arrows represent very weak interac- 
tion. Blue straight line: microwave alternating magnetic field pulse. 
Filled(empty) circle: initial(final) state of the electron. Setting the 
qubit state: i) excitation from the ground state to the ^Ai ex- 
cited state ii) spin-flip to the appropriate states iii) setting either 
Ms = 3/2 or i\fs ~ —3/2 state by tt pulse. Readout the qubit state: 
i) TT pulse to go back to Ms — ±1/2 states ii) spin-flip to the ^Ai 
state iii) radiative recombination to the ground state. We show 
the fine structure of states in the absence of external magnetic 
field (B=0) (2nd column) and at B > with typical EPR conditions 
in the next columns. D ~ 1685 MHz fine structure constant was 
determined by EPR f2lll . The figure does not show the true scale for 
the sake of clarity. 
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